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Pollen morphology was studied in 23 Iranian species of Cirsium (Asteraceae ) 
representing five sections, Cirsium, Cephalonoplos, Echenais, Epitrachys and 
Pseudoepitrachys by light and scanning electron microscopy. Among them five species 
are endemic to Iran and their pollen data are new to science. Pollen grains were radially 
symmetrical, isopolar, trizonocolporate and with echinate sculpturing. The tectum is 
microperforate, microreticulate or microgranulate in sculpture. The shape of pollen 
grains is suboblate to oblate-spheroidal in the equatorial view and is circular to triangular 
in polar view. The spines are conical with broad base and pointed or rounded tip. The 
number of spine rows between colpi varies from 4 to 5 among the species. ANOVA test 
performed for pollen characteristics showed no significant difference except for equatorial 
axis length, indicating that pollen data may not be used for delimitation or discrimination 
of those sections. PCA analysis showed the importance of echinae characteristics along 
with equatorial axis length, distance between two colpi and polar axis length in species 
delimitation. Clustering methods and the ordination plot of the PCA analysis support the 
affinity in the sect. Pseudoepitrachys and the distinction of the sect. Cephalonoplos from 
other sections, but did not support the grouping of the remaining sections from pollen 
morphology. 

Key words: Cirsium, cluster analysis, Iran, PCA analysis, pollen morphology. 


The genus Cirsium Mill, is one of the largest 
genera in the family Asteraceae, consisting of 
over 250 species of spiny perennial, biennial 
or rarely annual plants, and is distributed in the 
northern hemisphere, spanning the subtropical 
to boreal latitudes mainly distributed in Europe, 
North Africa, Asia, and North and Central 
America (Segarra-Moragues et al. 2007, Yildiz 


and Dirmenci 2008). The highest species 
diversity is concentrated in the mountains of 
Southern Europe and in the Caucasus (Bures et 
al. 2004). 

The genus is attributed to tribe Cardueae 
subtribe Carduinae in Asteraceae (Susanna 
et al. 2006). About 28 Cirsium species have 
been reported from Iran in the Flora Iranica 
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Table 1. Cirsium species examined in this study 


Section/Species 

Sect. Psudoepitrachys 

C. spectabile DC. 

C. conges turn Fitch. & C. A. Mey 
C. pyramidale Bomm. 

Sect. Epitrachys 
C. bracteosum DC. 

C. bornmulleri Sint, ex Bomm 
C. strigosum (M. Bieb.) M. Bieb. 
C. ciliatum (Murrey) Moench 
C. osseticum (Adams) Petrak 

C. aduncum Fitch. & C. A. Mey 
C. haussknechtii Boiss. 

C. lappaceum M. Bieb. 

C. turkestanicum (Regel) Petrak 
C. vulgare (Savi) Ten. 

Sect. Echenais 

C. echinus (M. Bieb.) Hand.-Mazz. 
Sect. Cirsium 

C. obvallatum (M. Bieb.) M. Bieb. 
C. libanothicum DC. 

C. cereticum (Lam.) d'Urv 
C. alatum (S. G. Gmelin) Bobrov 
C. elodes M. Bieb. 

C. canum (L.) All. 

C. hygrophilum Boiss. 

C. rhizocephalum C. A. Mey. 

Sect. Cephalonoplos 
C. arvense (L.) Scop. 


Locality and voucher specimen 


Kerman: Ghgher (Baft), Ghonbadan, Mirtadzedini 32098 (TUH) 

Tehran : Shahrestanak, Khor village, 2500 m, Nouroozi & Fathollahi 8600155 (HSBU) 
Kerman: Baft, Hankuye, 2300 m, Mirtadzedini 23104 (TUH) 

Kerman: Ghogher, Cheshme sabz, Mirtadzedini 29488 (TUH)(Endemic) 

Khorasan: 47 km from Shirvan to Sarani,1740 m, Ghahreman & Attar 27658 (TUH) 
Khorassan: Mashhad, 75 km to Kalat-e Naderi, Taherabad, 975 m, Ghahreman & Attar 27679 (TUH) 
Azarbaijan: Salmas to Khoy, 1400 m, Nouroozi 8600160 (HSBU) 

Azarbaijan: Arasbaran protected area, inter veighan and mahmoodabad, 1700 m, Ghahreman & 
Attar 17588 (TUH) 

Azarbaijan: Meshkinshahr, Mazrae Jahan,1169 m, Nouroozi & Fathollahi 8600170 (HSBU) 
Azarbaijan: Tabriz, Kandovan, Hilevar village,1600 m, Nouroozi 8600180 (HSBU) 

Azarbaijan: Meshkinshahr, Lahrood, Shabil, 2850 m, Ghahreman & Mozzafarian 70123 (TUH) 
Khorasan: Between Shirvan to Sarani, zivar village, Ghahreman & Attar 27661 (TUH) 

Tehran: Road of Karaj to Chaloos, Asara, 2050 m, Nouroozi 8600218 (HSBU) 

Azarbaijan: Marand, Yam village, Mishodagh, 2200 m, Nouroozi 8600185 (HSBU) 

Azarbaijan: Meshkinshahr, Lahrood, Shabil, 2850 m, Nouroozi 8600175 (HSBU) 

Azarbaijan: Marand, Yam village, Mishodagh, 1900 m, Nouroozi 8600190 (HSBU) 
Azarbaijan: Road of zabol to zahedan, Khaje mountains, 1300 m, Ghahreman (TUH) 
Azarbaijan: Jolfa, road of Kaimakidagh, 736 m, Nouroozi 8600150 (HSBU) 

Azarbaijan: 35 km from Kaleybar to Ahar, 1700 m, Ghahreman & Mozzaffarian 17534 (TUH) 
Azarbaijan: Khoy, 1500 m, Nouroozi 8600120 (HSBU) 

Tehran: Road of Haraz, lasem village,1900 m, Nouroozi & Fathollahi 8600135 (HSBU) 
Azarbaijan: Meshkinshahr, Lahrood, Shabil, 2850 m, Ghahreman & Mozzafarian 70121 (TUH) 

Azarbaijan: Bostanabad,1600 m, Nouroozi 8600120 (HSBU) 


(Rechinger 1979) and have been classified 
into five sections and several subsections and 
series in the Flora. The genus is considered 
to be taxonomically complex due to the 
variability and intergradations of diagnostic 
characters among taxa. The complex nature of 
the genus may be partly due to a high degree of 
hybridization occurring among different species 
and also due to highly variable chromosome 
numbers seemingly interacting with different 
environmental conditions. Several hybrids have 
been described from the Caucasus and adjoining 
regions of Asia Minor (Ownbey 1951, Moore 


and Frankton 1962, Charadze 1998, Bures et al. 
2004). 

Literature available on Asteraceae pollen is 
comprehensive, but publications have been very 
limited in Iran (Jafari and Ghanbarian 2007). 
Therefore, the present study analysed pollen 
morphology in 23 Cirsium species growing in 
Iran for the first time. 

Materials and Methods 

Pollen morphology was studied by light 
microscopy (LM) and scanning electron 
microscopy (SEM) in 23 Cirsium species, 
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Fig. 1. Representative LM micrographs of pollen grains in Cirsium species studied. A. C. spectabile 
(equatorial view). B. C. spectabile (polar view). C. C. pyramidale (equatorialview). D. C. pyramidale 
(polar view). E. C. congestion (polar view). F. C. congestum(e quatorial view). G. C . bracteosum 
(equatorial view). H. C. bracteosum (polar view). I. C.burnmeulleri (equatorial view). J. C. 
burnmeulleri (polar view). K. C. strigosum (equatorial view). L. C. strigosum (polar view). Scale bar 
= 10 pm. 
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representing five sections, growing in Iran (Table 
1). The voucher specimens have been deposited 
in the Shahid Beheshti University Herbarium 
(HSBU) and Central Herbarium of Tehran 
University (TUH). 

The pollen samples were obtained mostly 
from freshly collected specimens and also from 
dry herbarium materials. Fully matured anthers 
were removed from the specimens and prepared 
by the standard acetolysis method (Erdtman 
1952, 1969), and then they were mounted in 
glycerin jelly and sealed with paraffin wax for 
light microscopy. Morphological studies were 
performed on minimum 20 pollen grains for 
each species by an Olympus light microscope 
Model BH-2. 

For scanning electron microscopy, the 
acetolysed pollen grains were attached to 
aluminum stubs with double sided cellophane 
tape, air-dried at room temperature and coated 
with gold. The specimens were examined with a 
Philips XL 30 or a Leica LEO 440i at 15 kV and 
22 kV and photographed. 

UTHSCSA Image Tool ver. 3 software was 
used for pollen measurements, and then data 
obtained were scored. For multivariate analyses 
the mean of quantitative characters were used, 
while qualitative characters were coded as 
binary/multistate characters. Standardized 
data (mean = 0, variance = 1) were used for 
multivariate statistical analyses. The average 
taxonomic distance and squared Euclidean 
distance were used as dissimilarity coefficient in 
cluster analysis of data (Podani 2000). 

Palynological characteristics studied 
included, polar axis, equatorial axis, ratio of 
polar axis/equatorial axis (P/E), spine length, 
exine thickness, spine base width , number 
of spine row between colpi, mesocolpium 
length, mean distance between two spines, 
number of spines per 100 pm 2 , pollen shape, 
exine ornamentation and spine tip shape. The 
terminology is according to Punt et al. (2007). 

Analysis of variance (ANOVA) was 
performed among different sections studied for 


quantitative palynological characters noticed, 
to show if these data may be used in section 
delimitation. Principal component analysis 
(PCA) was performed among the species studied 
to determine palynological characteristics useful 
for separating the species. In order to group 
the species, cluster analysis using UPGMA 
(Unweighted Paired Group with Arithmetic 
Average) and NJ (Neighbor-joining) methods 
and PCA ordination plot were performed using 
euclidean and taxonomic distances calculated 
among the species (Podani 2000). ANOVA and 
least significant differences (LSD) tests were 
conducted using SPSS ver. 9 (1998), while 
clustering and ordination plot analyses using 
NTSYS ver. 2 (1998). 

Results 

Results of palynological study are presented 
in Table 2 and Figures 1-5. All species studied, 
showed trizonocolporate (Fig. 1), isopolar and 
radially symmetrical pollen grains; with echinate 
sculpturing. Pollen outline is suboblate or oblate- 
spheroidal in equatorial view and circular to 
semi angular in polar view. The mean of polar 
axis ranges from 61.93 pm in C. arvense (Fig. 
2D) to 79.85 pm in C. osseticum ; the equatorial 
axis ranges from 68.20 pm in C. aduncum 
(Fig. 2A) to 87.99 pm in C. obvallatum. The 
mean width of mesocolpium ranges from 45.41 
pm in C. alatum (Fig. 2B) to 71.75 pm in C. 
libanoticum (Fig. 2L). The exine ornamentation 
is echinate with microperforated tectum in most 
of the species studied (Fig. 3B-F, H) and with 
microgranulated in C. aduncum (Fig. 3A) and 
C. haussknechtii or microreticulated tectum 
in C. alatum (Fig. 3G), C. cereticum and C. 
obvallatum. The spines are conical with broad 
base and rounded (Fig. 3A-D) or rounded tip 
(Fig. 3E-H). The number of spine rows between 
two colpi varies between 4-5, resulting into 
about 12 to 15 spines in the equatorial plane. The 
distance between two spines ranges from 15.89 
pm in C. aduncum (Fig. 2A) to 24.03 pm in C. 
cereticum. The mean number of spines per 100 
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Fig. 2. Representative SEM micrograph of pollen grains in Cirsium species studied. A. C. aduncum. B. 
C. alatum. C. C.bracteosum. D. C. arvense. E. C. spectabile. F. C. burnmeulleri. G. C. ciliatum. H. C. 
rhizocephalum. I. C. pyramidale. J. C. lappaceum. K. C. hygrophilum. L. C. libanoticum. 
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Fig. 3. Representative exine ornamentation in Cirsium species studied. A. C. aduncum. B. C. 
spectabile. C. C. arvense. D. C. bracteosum. E. C. burnmeulleri. F. C. pyrcimidale. G. C. alatum. 
H. C. osseticum. 













Table 2. Pollen morphological parameters in Cirsium species examined. Data of statistical columns consists of the mean ± standard deviation in top line and rang in bottom line ximumon bottom line . 
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Fig. 4. UPGMA dendrogram of palynological data in Cirsium species studied. 


pm 2 varies from 2 in C. libanoticum (Fig. 2L) to 
11 in C. aduncum (Fig. 2A). The mean of spine 
length varies from 6.50 pm in C. aduncum (Figs. 
2A, 3A) to 10.00 pm in C. spectabde (Figs. 2E, 
3B). The mean width of spine base varies from 
9.48 pm in C. aduncum (Fig. 3A) to 16.13 pm in 
C. canum. The mean of exine thickness ranges 
from 8.56 pm in C. arvense to 11.59 pm in C. 
obvallatum. 

ANOVA test performed for pollen 
characteristics among five sections of Iranian 
Cirsium studied showed significant difference 
only for equatorial axis length (F = 3.14, p 
< 0.05). However, PCA analysis performed 
on these characters revealed that the first 3 
components comprise about 65% of the total 
variation and in the first component (38% of 
variation), characters like equatorial axis length, 
polar axis length, echinae length, distance 
between two echinae, echinae base width and 
distance between two colpi have the highest 
positive correlation (r > 0.70). 


UPGMA, NJ and WARD methods of 
clustering of pollen data produced similar 
results- with a higher cophenetic correlation (r = 
0.95), which was also supported with respect to 
some species by ordination plot obtained (Figs. 
4,5). 

Considering the species affinity represented 
in the UPGMA dendrogram (Fig. 4), in sect. 
Psudoepitrachys C. congestion, C. pyramidale 
and C. spectabile formed a subcluster, which 
shows close affinity to each other. In sect. 
Epitrachys, C. bracteosum and C. turkestanicum 
were placed close to each other. Cirsium 
haussknechtii then joined these species with 
some distance. Cirsium ciliatum, on the other 
hand, showed affinity with C. strigosum, and 
C. osseticum joining these two species at 
some distance. However, C. aduncum and C. 
iappaceum were placed in different subclusters, 
respectively. In sect. Cirsium, C. elodes and C. 
rhizocephalum formed a subcluster showing 
some degree of similarity, and C. canum joined 
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Dim-1 


^ Sect. Epitrachys □ Sect. Cirsium ^ Sect. Cephalonoplos 

A Sect. Psudoepitrachys ^ Sect. Echenais 


Fig. 5. PCA plot of palynological data in Cirsium species studied. 


these species with some distance. Cirsium 
libanoticum and C. hygrophilum as well as C. 
cereticum and C. obvallatum, however, showed 
affinity to each other but they were placed 
in different subclusters. Cirsium alatum was 
placed in the most separated subcluster with C. 
aduncum of the sect. Epitrachys. 

PCA analysis showed that the first three 
components (coordinates) comprise about 65% 
of total variance in which the first, second and 
third components explained about 38%, 15% 
and 12% of the total variation, respectively. The 
ordination plot on first two component axes 
showed that plots of the sect. Psedepitrachys and 
the sect. Cephalonoplos are separated from those 
of other sections along Dim-2 (Fig. 5). 

Discussion 

The species investigated in the present study 
are from all sections of Iranian Cirsium. Pollen 
grains observed are usually radially symmetrical, 
isopolar, trizonocolporate, oblate spheroidal 


to spheroidal with echinate sculpturing, which 
is almost in agreement with previous reports 
on the genus Cirsium, although the shape of 
pollen grains was slightly differently described 
as suboblate or oblate-spheroidal (Halbritter 
et al. 2005), spheroidal to oblate-spheroidal or 
spheroidal to prolate-spheroidal (Meo 2005), 
spheroidal (Hesse et al. 2009), spherical (Jafari 
and Ghanbarian 2007), subprolate or prolate 
spheroidal (Yildiz et al. 2011) in equatorial view 
and is circular (Ylldiz et al. 2011) or semiangular 
or semiangular to circular in polar view (Meo 
2005). 

According to Flora Iranica (Rechinger 1979) 
the genus Cirsium has five sections including 
Cirsium, Cephalonoplos, Echenais, Epitrachys, 
and Pseudoepitrachys. In the present study, the 
species attributed to sectt. Pseudoepitrachys and 
Cephalonoplos are all suboblate, while those 
attributed to sections Epitrachys and Cirsium are 
oblate-spheroidal. All the species have echinate 
sculpturing, and the exine ornamentation 
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between spines is microperforated in all 
species of sect. Pseudoepitrachys, Echenais 
and Cephalonoplos. The species of sect. 
Epitrachys have also microperforated exine 
excluding C. aduncum and C. haussknechtii 
that have microgranulated exine. Cirsium 
aduncum and C. haussknechtii are separated 
and attributed to subsect. Microcephala in the 
sect. Epitrachys due to the presence of small 
capitula (Rechinger 1979). The species of sect. 
Cirsium have microperforated exine excluding 
C. cereticum, C. alatum and C. obvallatum that 
have microreticulated exine. 

The ANOVA result indicated that five 
sections of Iranian Cirsium do not differ 
significantly in their pollen characteristics and 
such data may not be used in discriminating 
different sections. However, the PCA result 
which represented the first three factors 
comprising about 65% of all of the variation, 
showed that characteristics such as polar axis 
length, equatorial axis length, echine length, 
distance between two colpi, distance between 
two echinae and the echine base length are the 
most variable pollen characteristics in the first 
factor with about 38.41% of the total variation 
and these characteristics may be used for species 
deli mi tation. 

Clustering and ordination plot of the species 
based on pollen data also indicated that almost 
sectional delimitation is not possible by such 
characteristics because the species from different 
sections have been placed close to each other 
and are intermingled. However, in some places 
good separation has occurred, for example the 
species of C. congestum, C. pyramidale and 
C. spectabile from the sect. Psedepitrachys 
almost form a separate group in the ordination 
plot (Figs. 4, 5). Moreover we find that in the 
sect. Epitrachys, C. strigosum, C. ciliatum and 
C. osseticum show close affinity supporting the 
treatment in Flora Iranica with respect to these 
species (Rechinger 1979). In sect. Cirsium the 
affinity of C. cereticum and C. obvallatum and 
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that of C. elodes and C. rhizocephalum, both 
suggested in Flora Iranica (Rechinger 1979), 
were also supported by their similarity in some 
degree from pollen characteristics. 

Our previous cytotaxonomic study on 
Iranian Cirsium species (Nouroozi et al. 2010) 
revealed the affinity between C. congesum 
and C. spectabile from sect. Pseudoepitrachys 
and that among C. alatum, C. hygrophilum, 
C. libanoticum and C. obvallatum from sect. 
Cirsium. Both affinities are almost supported by 
the present palynological data. 
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